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We present a theory of negative Coulomb drag in capacitively coupled quantum wires based on 
the commensurability of the electron density and the long-range nature of the Coulomb interaction 
in the Tomonaga-Luttinger liquid. The commensurability introduces a notion of doped particles 
and holes. We point out that the long-range interaction allows a particle-hole pairing over the wires 
and that the particle-hole pairing brings about the positive drag of holes, that is, the negative drag. 
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Understanding many-body interaction effects is a long¬ 
standing objective of condensed matter physics as typi¬ 
fied by the Fermi liquid theory [I]. The Fermi liquid 
theory however fails in strongly correlated systems, es¬ 
pecially in one-dimensional (ID) systems [5j. Since any 
particle in ID space cannot overtake others in front of 
it, the restricted ID geometry provokes characteristic 
transport phenomena such as the quantization of conduc¬ 
tance 0111 and diffusive transports compatible with in- 
tegrability [Mi. In ID systems, instead of the Fermi liq¬ 
uid, the collective motion of the particle-hole pair in the 
so-called Tomonaga-Luttinger liquid (TLL) is the most 
elementary excitation mm- 

The TLL has a great advantage of incorporating in¬ 
teractions into a single parameter K known as the TLL 
parameter. It is non-interacting for K = 1, repulsive for 
K < 1 and attractive for IF > 1 [2j. Furthermore K 
is controllable with external parameters. In the field of 
quantum magnetism, the experimental controllability of 
K with the magnetic field is fully exploited to simulate 
itinerant boson systems mm- Clearly this idea of exter¬ 
nal control of interaction strength of TLL fits well with 
mesoscopic physics. Nevertheless it is not fully empha¬ 
sized thus far. 

One of representative mesoscopic systems yielding the 
TLL is the quantum wire. With the aid of the long- 
range Coulomb interaction, the TLL parameter of the 
quantum wire is expected to be small. Such a strongly 
repulsive TLL, or the ID Wigner crystal puna , in quan¬ 
tum wires is a highly nontrivial electron state thanks to 
the interplay of the long-range interaction and the non- 
Fermi-liquid nature. 

Recent progresses in fabricating quantum wires have 
made it possible to address experimentally a remarkable 
transport property of TLL, that is, negative Coulomb 
drag HME5]. Coulomb drag in general is an induction 
phenomenon of the electric current by another capaci¬ 
tively coupled current (Fig. [TJ) , purely originating from 
the long-range Coulomb interaction. Coulomb drag of 
quantum wires has received intensive theoretical inter¬ 
ests m - f^l for decades. These theories predicted the 


positive drag of parallel drive and drag currents. How¬ 
ever, in sharp contrast to the theoretical predictions, the 
negative drag of antiparallel drive and drag currents was 
observed experimentally [T3]. The negative drag was 
“unexpected” m in this sense and is still elusive because 
of lack of theoretical explanation. The authors of Ref. [13; 
indicated that the Wigner crystal formed on the drag 
wire would be responsible for the negative drag. Later 
Refs. 12] and [13 reported that an increase of the voltage 
induces alternately the positive and negative drags. The 
Wigner crystallization on the single wire is not likely to 
explain this re-entrant negative drag [14] . Therefore, the 
mechanism of the negative drag remains unclear despite 
the existence of the fascinating experimental results. 

In this Paper, we propose a simple theory that explains 
the negative drag clearly. Its mechanism is summarized 
as follows. When the particle density of the drive wire is 
commensurate with the hole density of the drag wire, the 
particle-hole pairing over the wires occurs, leading to the 
positive drag of the hole current , that is, the negative drag 
of the electric current as shown in Fig. [3] (b). We point 
out that the strong repulsion of the TLL is necessary to 
induce the negative drag. We also clarify the crucial role 
of the controllability of interactin strength in the negative 
drag. 

Let us consider quantum wires of the length L and the 
width d w separated by the distance d (Fig. [l]) and apply 
the external voltage Vdrive to the drive wire in order to 
induce the drive current /drive- Our purpose is to see the 
sign of the ratio Idrag//drive- If the interwire interaction 
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Figure 1. Coupled quantum wires (the shaded area). 
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was absent, each wire would have a TLL 


U n = 


Hu n 
2i r 


dx 


K n {d x e n ) 2 + — (d x fa) 2 

An 


( 1 ) 


Here u n is the velocity of the TLL and K n (< 1) is the 
TLL paramter. We label the drive and drag wires by 
n = 1 and n = 2. fa and 9 n satisfy the commutation 
relation [fa(x),d x 9 m (y)\ = iTi8 n ^ m d(x-y), leading to the 
equation of motion dtfa = u n K n d x fa. The Hamiltonian 
(|T|) describes the bosonic excitation with the linear dis¬ 
persion, that is, the TLL. Since the TLL describes the 
particle-hole excitation, the bosonic field <p n is related 
to the fluctuations of the electron density p n (x) from its 
average p n in the nth wire [7 , 


d (f) ^ 

Pn{x) = pn - —~ + 2p n V' cos(2p7rp n x + 2 pfa). (2) 

7r z ' 

p -1 


The interwire interaction induces two effects: the 
renormalization of the TLL parameters K n and a lock¬ 
ing of bosonic fields over the wires. The latter needs 
a careful consideration. Let us consider the interwire 
Coulomb interaction f drdxV±(r)pi(x)p 2 (r + x ). Al¬ 
though the Coulomb interaction is originally long-ranged, 
it is replaceable to the effective short-range interaction 
gpi(x)p 2 (x) in the low-energy limit. The long-rang na¬ 
ture of the interaction renormalizes the strength g of the 
effective interaction. We will come back to this point 
later. The short-range intearction gpi(x)p 2 (x) generates 
an interaction Agpifa cos(2pi7rpi2;+2pi</>i) cos(2pi7rp2:E-|- 
2 p 2 ^ 2 ), that is, 


y; 2gpiP2COs[2ir(p 1 p 1 +vp2P2)x+2(Pi<t , i+vP2(f>2)], (3) 

v=± 


where p\ and P 2 are integers. In the case of a single 
quantum wire, these cosine interactions are irrelevant ex¬ 
cept for the Mott insulating case with repulsive interac¬ 
tion, when p = 1/ao, where a o is the lattice spacing of 
the wires [2]. However, in our model of coupled quan¬ 
tum wires, p± and p 2 are independently tunable param¬ 
eters. This is important for the negative drag as we dis¬ 
cuss below. In the following we consider the case with 
Pi = P 2 = P- It is easy to extend the results for general 
cases with pi ^ p 2 . 

Including these effects, we consider a model with a 
Hamiltonian 



K„(d x O v ) 2 


1 

Kl 


(d x fa) 2 



( 4 ) 


where fa = (fa ± fa)/\/2, 9± = (6i ± 6» 2 )/v / 2, 


VP = 2g PlP 2 / dx cos(VpTt Pl ,x + pV8(f> v ), 
Jo 


(peN), 

( 5 ) 


and p± = (pi ± p 2 )/ 2. 

When p v is commensurate, that is, 


2pa 0 p u e Z, (6) 

the interaction © survives the spatial integration. Oth¬ 
erwise the interaction |5| is negligible because of the de¬ 
structive oscillation . In the commensurate case 

([b]), the cosine interaction ([5]) causes an energy gap for 
t fa and fixes fa to a constant value. This is called the 
locking effect. 

Electric currents of the symmetric and antisymmetric 
sectors are given by J± = —e(d t <f>±)/ 7r, where — e < 0 is 
the electron charge. Similarly one can define the current 
J n = —e(dtfa)/ 7r on the nth wire. As well as the in¬ 
commensurate density, the incommensurate current con¬ 
flicts with the locking effect because the current adds a 
temporal oscillation e lJ v nt / e to V£. Thus, if V? eventu¬ 
ally locks the fa field over the wires to a constant, then 
= 0 follows [24]. Note that two currents J± do not 
vanish simultaneously. If J+ = J_ = 0, the circuit has 
no current, which is different from the situation that we 
consider here. The drag and drive currents are given by 


-‘drag 


J+ - J 


^ drive 


J -f H - J— 


( 7 ) 


J-f = 0 leads to the negative drag Idrag/fdrive = — 1 and 
J_ = 0 leads to the positive drag /drag/^drive = 1. This 
is the essence for the negative drag. If we take p\ ^ p 2 
in Eq. <©> the locking of fa oc (fa ± (p 2 /pi)fa) leads to 

Alrag/Llrive = Vl /P i ■ 

Physical picture becomes clear when we take the 
particle-hole description. We assume pi > p 2 and 


V' 

P+ = o-• 

2 pa 0 


( 8 ) 


p and p' £ N are coprime. Compared to the average 
filling ([8]), the drive (drag) wire has more (less) electrons. 


(a) positive drag 


E 



(b) negative drag 
E 
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Figure 2. (a) The situation (101 when the positive drag oc¬ 
curs. Efti is the Fermi level of the ?ith wire. The drive and 
drag wires have the equal particle densities, (b) The situa¬ 
tion when the negative drag occurs. The density of doped 
particles in the drive wire is equal to that of doped holes in 
the drag wire. 
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Thus one may regard the drive and drag wires as particle- 
doped and hole-doped conductors respectively. In this 
particle-hole description, the condition ([6]) for v = + 
becomes 


p\ = p\, (9) 

where = p n — p'/2pao and p 1 / = — pP are densities 
of doped particles and holes in the nth wire. Thus one 
can see that the negative drag occurs when the density 
of doped particles and holes are balanced. Similarly the 
positive drag originates from a commensurate condition 
P- = 0, that is, 


P P 

Pi = Pi¬ 


rn 


We visualized the conditions © and ( [lOt Fig. [2] When 
the condition (10) holds, the strong interwire interaction 


locks relative positions of particles over the wires to avoid 
costing the Coulomb potential energy VZ, resulting in the 
positive drag [Fig. [3] (a)]. On the other hand, when the 
condition © holds, the strong attraction works between 
particles on one wire and holes on the other wire, form¬ 
ing particle-hole pairs over the wires. This particle-hole 
pairing allows for the positive drag of the hole current, 
that is, the negative drag [Fig. 1(b)]- 

As shown above, the simple model Q explains the pos¬ 
itive and negative drags on equal footing. The originality 
of the present model is in the inclusion of the cosine inter¬ 
action V? of the symmetric sector. Except this point, the 
model Q is identical to the models commonly employed 
in theoretical studies of Coulomb drag urmaiza. The 
preceding theories implicitly assumed the incommensu¬ 
rate p + and ignored V+ [TBH22| . This assumption is cer¬ 
tainly natural for ID conductors, but not always true. 
There is no principal reason to rule out the possibility 
([8]). In the following we discuss the situation where the 
condition Q holds. 

Since quantum wires are usually far below the unity 
filling (i.e. p n <C 1/ao), large p is requisite for satisfying 
Eq. ([8]). In this case, the relevance of the cosine inter¬ 
action is reduced. This is because the cosine interaction 
V? has a scaling dimension d = 2 p 2 K+, and d < 2 is 


(a) positive drag (b) negative drag 



Figure 3. Cartoons of (a) the positive and (b) the negative 
drags, (a) The interwire repulsion of particles results in the 
global shift of the particles, that is, the positive drag, (b) 
The negative drag is the positive drag of holes due to the 
particle-hole pairing. 


required in order to lock <j> + . Larger p makes it more 
difficult to achieve d < 2. In contrast the positive drag 
does not require the large p because the condition ([h]) is 
easily satisfied for p\ = p 2 . This point makes the pos¬ 
itive and negative drags unequal. The strong repulsion 
K + < 1/p 2 <d 1 is necessary for the negative drag. 

Let us show that the long-range intrawire and inter¬ 
wire Coulomb interactions are the keys for realization of 
the condition K + < 1/p 2 -C 1. The intrawire Coulomb 
interaction f dxdrVjj(r)p n (x)p n (x + 7') renormalizes u n 
and K n as follows. It generates the kinetic term Vj| (q = 
0){d x (f> n ) 2 /n 2 , where V\\(q) is the Fourier transform of 
Vj|(r). When the electron density is low enough, the 
Coulomb interaction is unscreened and divergent at q —»■ 
0. Such a divergence at long distance is crucial to the 
drag current flowing through the whole wire. In order 
to include the effect of the cut-off of the divergence due 
to the finite length of the wires, we model the intrawire 
interaction as Vy(r) = e 2 /47T£yV 2 + d 2 [25] ; where £ is 
the dielectric constant of the wires. Because of the di¬ 
mensionality, Vy(g) ~ (e 2 /27T£) log(l/|g|d w ) exhibits the 
logarithmic divergence IH1251. The finite length L cuts 
off the divergence as V\\(q —► 0) ~ (e 2 /27re) log(L/d w ), 
leading to [24] 


Kr, ~ K = 


7 T 2 h,£VF 


log (L/d* 


- 1/2 


( 11 ) 


and u n = u = vf/K. Here vf is the Fermi velocity 
of the free electron. Since the Coulomb interaction in 
ID diverges at both long and short distances, the TLL 


parameter (111 depends on both L and d w . The relation 


(111 shows that we can control the strength of repulsion 


of the TLL by changing geometrical parameters L and 
d w of the circuit. One can find a similar argument in 
Refs. [26] and [27. 

The interwire interaction is similarly given by V± (q —> 
0) ~ (e 2 /27T£') log (L/d), where s' is the dielectric con¬ 
stant of an insulating medium in between the wires. The 
effective q independence of V±(q —> 0) allows us to re¬ 
place it to an effective short-range interaction V±(r) — 
V±(q —> 0)<5(r) in the low-energy limit q 0, which gives 
the strength g = V± (q —> 0) of the cosine interaction (J5|. 

Thus both of the intrawire and the interwire interac¬ 
tions renormalize the TLL parameter as [23] 


I<± = 


j^[log(L/i.)±^ ; log(L/d) 


1 - 1/2 


( 12 ) 


and the velocity to u± = vp/K±. In general K + < K_ 
holds. The TLL parameters (12) are controllable with the 
geometrical parameters of the circuit. When the wires 
are long (L d w ) and separated by a medium with 
e' /?> £, The TLL parameter K + ~ [log(L/d w )] -1 / 2 is 
basically determined only from the intrawire Coulomb 
intearction. Since K + ~ [log(L/d w )]” 1//2 can become 
arbitrally small, the condition K + < 1/p 2 <C 1 for the 
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negative drag can be satisfied with the aid of the long- 
range intrawire Coulomb interaction. 

Let us discuss major factors that disturb the negative 
drag. They are the temperature and the incommensu¬ 
rability. First we estimate the temperature effect. Un¬ 
der the condition ([ 8 ]), the symmetric mode acquires an 
energy gap for <j> + . We can easily see that an expan¬ 
sion cos(p\/ 8 </>-|_) ~ 1 — 4 p 2 (4>+ — Q\) 2 around a bottom 
tf> + = Q+ of the cosine V+ generates the quadratic mass 
term [24] . However the locking is weakened by instan- 
tons that represents the tunneling of neighboring locking 
values of Q° + —► Q° + ± Tr/py/2. Therefore, in order to 
have Coulomb drag, the excitation of instantons should 
be suppressed nans]. Fortunately the exact excitation 
spectrum of the instanton is available |28l I29| . which is 
composed of a soliton and an antisoliton. Their excita¬ 
tion gap, M, is exactly derived [24] [30]. We can suppress 
thermal excitations of the instanton at 


k B T < M 


2e 2 p+ 

pn 2 £ 


2 e 

e' 


1/2 


log(L/d w )log(L/d) 


(13) 


Then 4> + is well locked to the bottom Q+. As we dis¬ 


cussed, the negative drag occurs when log(L/d w ) 1 
and e/e' <C 1 are satisfied. Given log(L/d w ) s'/e 1, 
the gap M becomes large and then the temperature range 


(13) is wide enough for experimental realizations. 


Next we estimate robustness of the negative drag 
against the incommensurability. Let us displace p+ from 
the commensurate value We use = Apn\p + — p\_\ 
as a measure of the displacement. Nonzero adds 
the incommensurate oscillation e lx ^ p to the cosine in¬ 
teraction ([5]), which disturbs the coherence of Q+. Let 
us write the lowest-energy excitation gap of </>+ as u+/£ 
( 7 ^ M) [24]. £ gives the coherent interval of the locking 
of </>+. When £ < £ p , the oscillation is very slow over the 
length £ , and the locking by V+ persists. On the other 
hand, when £ > £ p , the oscillation dissolves the locking. 
Therefore the negative drag lasts for £~ 
is El, 


\p+~P° + \ 


P+ 


< A = n 


2 e log (L/d) 
_£' log (L/d w )_ 


< r\ 

that 

1/2 

(14) 


for e'/e > 1 . 

Furthermore, in order to realize the negative drag, the 
positive drag should be suppressed. Since the positive 
drag is stabilized for 

( 15 ) 

P+ 


for e'/e 1 , the electron densities must break the in¬ 
equality (15). Thus the negative drag occurs when pi 
and p 2 satisfy | p+ - p+\ < Ap+ < \pi - p 2 \- 

Let us compare our theory with the existing experi¬ 
ments. According to and (fl2|, we can prepare K± 


as small as we wish by taking large enough L/d w . In 
fact an experiment shows A"_ ~ 0.08 ± 0.02 [15]. Such 
a small TLL parameter is impossible without the long- 
range repulsion. For instance, the ID Hubbard model 
only with the on-site repulsion has the TLL parameter 
larger than 1/2 for any filling and parameters [31]. The 
extraordinary small TLL parameter observed experimen¬ 
tally indicates that the long-range nature of the Coulomb 
interaction does exist in quantum wires. 

In general the ./V-component TLL on the quantum wire 
has a quantized conductance G = Ne 2 /h [3]. When the 
spin is included, G = N'(2e 2 /h) with N = 2TV'. In our 
model, the interaction V/j locks a half degree of freedom 
in </>2 = — (f>-)/V 2, resulting in a fractionalization 

G = (1/2 )e 2 /h, similarly to Ref. EU Inclusion of the 
spin degree of freedom doubles it to G = e 2 /h, which is 
basically consistent with experiments m- However the 
height is not well quantized as e 2 /h for some situations: 
Fig. 3a of Ref. [14] shows that the height of conductance 
plateaus are reduced. Such reduction is attributed to, 
for example, the impurity [5] and the spin degree of free¬ 
dom [271 and irrelevant in the emergence of the negative 
drag. Our theory explains at least the negative drag on 
the lowest conductance plateau observed in Fig. 4a of 
Ref. 1111 Those on higher plateaus can be explained after 
a straightforward extension. When intrawire intercom¬ 
ponent interactions are negligible, which is usually true 
in the low-energy limit, we can easily extend our model to 
the Wcomponent TLL on higher conductance plateaus. 

Figure [2] implies that the positive and negative drag 
occurs alternately as pi and p 2 are increased indepen¬ 
dently. Besides, Coulomb drag becomes less promi¬ 
nent in larger p + . Large enough p + easily satisfies 
maxll/U — p\\, \p\ — P 2 1} < Ap+, inducing the positive 
and negative drags simultaneously [Eqs. ( fl4] ) and ©], 
that is, no drag as a whole. Such qualitative dependences 
on the densities p\ and p 2 are consistent with the exper¬ 
iments mm by translating gate voltages to electron 
densities. 
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Supplemental Material for 
“Negative Coulomb Drag in Coupled 
Quantum Wires” 


LOCKING OF THE ZERO MODE 


Here we show that the locking of </>+ leads to ( d t (j)+) = 
0. To see this, we introduce a mode expansion [2j, 

<t>+(x) = Q+ + —^~ x 


T.fm no 


0+{x) = 0+ H- —^x 

Ju 


^sgn(-g)y^±^ (q q , +e lqx + °l + e ^ x ). 

(17) 


The first lines of Eqs. (16) and ( [Tt) ) represent the 
zero modes of 4>+{x) and 0+(x), where ((5+,0+) and 
(P + ,N + ) are canonical conjugate, that is, [Q + ,P+] = 
[0 + ,iV-|_] = i. The second lines represent the g / 0 
modes. In the absence of the cosine interaction, a q:+ is 
an annihilation operator of the TLL. Under the condi¬ 
tion (6) of the main text, the effective Hamiltonian of 
the symmetric mode becomes 


n + = 


hu + 

27r 


dx 


K + (d x 9+) 2 + — (cW>+) 5 
K + 


2gpiP2 / dx cos(pV8(f)+) 


(18) 


with 


g = v±(q 


° )= 2)rP l0S(L/d) - 


The mode expansions (16) and (17) lead to 

7rfi/u + 


n + = 


- Nl 

2 K+L + 


nHu+K + 
2 L 


(19) 


P+ +^2hu + \q\al + a q , + . 

q^O 

( 20 ) 


In the presence of the cosine interaction, the TLL ac¬ 
quires the excitation gap. When the temperature is much 
lower than the energy gap of the soliton, we can expand 
the cosine as cos(p-\/8<)> + ) ~ 1 — 4p 2 (</> + — Q+) 2 with a 

An I-bn TT r\ ill- /v n i n -»-» / 1 Q \ 


constant Q°, . Then the Hamiltonian (|18|) turns into 




TlhU-l 


-Ni 


■nhu+K. 


+ d2 


Pf + 


LMl 


2K + L~' + 2 L ' + 2-Khu + K_ 

+ ^x/O+l) 2 + (Ml /h) 2 b\ + b q ,+ . 

o 


(Q+ - Q +) 2 


(21) 
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We introduced the creation and annihilation operators, we obtain 
b q + and , through a Bogolioubov transformation, 


°q,+ 


cosh 9 q}+ — sinh 9 q ^ + 


— sinh 0 9 _ 


cos 9 q _ 


a q..+ 

2 ~q,+ 


( 22 ) 


with coth(20 g _) ~ 1 + 2(hu + q/Mi) 2 . Since there is no 
Bose-Einstein condensate, that is, ( 6 g ,+) = ( 6 g,+) = 0, 
the average ( d t <j)+) = u + K + (d x 9 + ) is fully determined 
from the zero mode, 


(M+) = ^±{ P+ ). 


(23) 


As one can see from Eq. (21), the zero mode P+ is subject 


to the Hamiltonian of a harmonic oscillator, 


•T/zero = nflU + N 2 , ^ hu + K + p2 

+ 2 K+L +2 L + 


LM ' (Q+-Q °+) 2 


2irhu-i r Kj r 


(24) 

Note that the zero mode N + is independent of both P + 


and Q + . Although Eq. (24) is, of course, the trivial con¬ 
sequence of the quadratic approximation of the cosine, it 
leads to an important results, 


(P+) = (( Q+ - Q°+)) = o. 


(25) 


(c>t</>+) = 0 immediately follows from Eq. (25). The rela¬ 
tion (25) holds only when the quadratic approximation 


of the cosine is valid. 
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where we used the features, T(l) = 1 and r(z)T(l — z) = 
7 r/sin (7 tz), of the Gamma function T(z). u±/K± and 
u±K± satisfy 
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and u±K± = uK = v F , resulting in 
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and u± = Vp/K±. In the end, when e'/e /g> 1, the gap 
of the soliton is approximated to 
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Then the lowest-energy excitation gap Mi is reduced to 

M\ = 2M sin(7T7/2) 
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Now it is straightforward to derive 


GAP OF THE SOLITON AND THE 
ANTISOLITON 
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The field theory with the Hamiltonian (18) is called 
as the sine-Gordon theory. Since the sine-Gordon theory 
is integrable, its various quantities are exactly available. 
The excitation gap of the soliton and the antisoliton at 
q = 0 is given by l-TD 
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For p 2 K + -C 1, the parameter 7 = l/[(p 2 K + ) 1 — 1] ~ 
p 2 K + is also small. Taking the lowest order of 7 ~ p 2 K +) 


The condition 1 < £ 1 in the main text can be rewrit¬ 
ten as 
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which becomes Eq. (14) in the main text. On the other 
hand, the corresponding condition for the positive drag 
is 
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